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Abstract 

In dense matter like the Supernovae and at high energy, neutri- 
nos have two points where the nonadiabatic transitions could occur. 
With the present values of the oscillation parameters in the standard 
three flavor scenario, two nonadiabatic transitions can be regarded 
as independent, i.e., two can be treated by the two flavor formalism 
separately. However, in the presence of new physics, it is not clear in 
general whether such independence holds. We examine this question 
by assuming hypothetical range of the neutrino energy and by varying 
the mixing angles and the mass squared differences in the standard 
case. We found the cases where Independent Crossing Approximation 
breaks down for some unrealistic range of Am^ and #13. We also 
discuss the criterion which gives such independence. 



1 



1 Introduction 



We now understand the solar neutrino problem pQ, i.e., we now know that 
the deficit of the electron neutrino flux from the Sun can be accounted for 
approximately as the two flavor neutrino oscillations in matter. From the 
precise experiments, we know the mixing angle of the neutrino oscillations 
is large. It is called Large Mixing Angle solution. Historically, the so-called 
Small Mixing Angle solution has also been discussed [2], because it gave 
good fit to the solar neutrino data in the past. In this case, the nonadiabatic 
transitions, which are jumping between the different neutrino energy levels 
in matter, would become important. On the other hand, we also understand 
now the atmospheric neutrino problem [Tj, i.e., the deficit of the muon neu- 
trino flux from the cosmic ray can be explained approximately as the two 
flavor neutrino oscillations. We also know the oscillation channel of the at- 
mospheric neutrino differ from the solar neutrino. These lead us to conclude 
that the neutrino oscillation is described by three flavor neutrino oscillations. 

In the case of the environments with high density, such as supernovae, it is 
important to consider the neutrino oscillations with the three flavor exactly, 
since there can be two energy level crossings. Furthermore, if nonadiabatic 
transitions occur at each level crossing at all, then we need to take them 
into account. In this case, there could be two nonadiabatic transitions [3]. 
The three flavor neutrino oscillations with two nonadiabatic transitions are 
very complicated in general. Although little is known about the exact so- 
lution of this very complicated problem, for instance, Ref. [I] assumes the 
Independent Crossing Approximation (ICA), which regards two nonadiabatic 
transitions as independent. In this approximation, we can regard each nona- 
diabatic transition locally as that between the two levels, and then we can 
get the total nonadiabatic transition just by multipling each probability of 
the nonadiabatic transition [I]: 

Pc -PLower ^ -^Higher 

The reason that ICA has been assumed is (i) because the two energy 
ranges which give a resonance do not overlap for the solar and atmospheric 
neutrino oscillation parameters and (ii) because the energy of the supernovae 
neutrino is so low and the solar neutrino mixing angle is so large that we 
would not have nonadiabatic transitions at both of the energy level crossings. 

However, if we have new physics beyond the standard model, then the 
standard matter effect may be modified [3 El [7] and it becomes unclear 
whether these reasonings always hold. The purpose of this paper is to exam- 
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ine whether ICA holds or not in the three flavor neutrino oscillations, and 
to discuss under which condition ICA breaks down. For simplicity we take 
only the three flavor case with the standard matter effect as an example and 
we will assume hypothetical neutrino energy and hypothetical values of the 
oscillation parameters to have nonadiabatic transitions at both energy level 
crossings. 

From numerical calculations, we found that ICA breaks down for some 
range of Am^ and #13. We found that ICA breaks down when the ratio 
Am^/Amlx is of order one and the mixing angle #13 is large. In the realistic 
range of the oscillation parameters, since Am^/Am^ is much larger than 
one and the mixing angle 6^3 is small, the two neutrino energy ranges for the 
resonance never overlap. This is the reason why the treatment in Ref. [I] is 
regarded as correct. However, it turns out that it is not the correct criterion 
because we find that two level crossings become independent even if the two 
resonances overlap. Instead of using the notion of the overlapping resonances, 
we will introduce the new parameters which gives the criterion of ICA. With 
this criterion of ICA, we interpret the breaking of ICA as the sign of a large 
contribution of the extra off-diagonal coefficients. Our numerical calculations 
confirm this analytic treatment very well. 

Although we discuss the simplest case only with the neutrino-electron 
interaction in this paper, we can apply our treatment to any other cases even 
if it is not known whether the nonadiabatic transitions are independent or 
not. For example, even if the effective ratio (Am^/Am^) in the presence of 
new physics is small, or the effective #13 is large, we can discuss independence 
of the nonadiabatic transitions. 

2 Neutrino oscillations with two flavor in mat- 
ter 

In this section, we review the standard treatment of the neutrino oscillations 
with two flavor. 

2.1 Adiabatic transitions 

The positive energy part of the Dirac equation for the flavor eigenstates 
v a (t) (a = e, fi) propagating at time t in matter is given by 




Hu a {t), 
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where H is given by 



h = u( ? !V' + f!I ), (2) 



OA/ \fl 

here A = Am 2 /2E, A = \f2G-pn e (R) is the extra potential in medium at a 
distance R = ct from the initial point, and U is the 2x2 MNS matrix [8]. 
can be diagonalized as 

H = exp(i<T 2 0)diag(£i, E 2 ) exp(-ia 2 0). (3) 

where a 2 is the Pauli matrix and the effective mixing angle 9 is given by 

• 2~n (Asill2#) 2 

qin 29 — - - (A) 

(Acos2fl- A) 2 + (Asin2#) 2 ' y> 

sin 2 29 becomes the maximum at the point A(R res ) = A cos 29, which is called 
the Resonance point. 

In the adiabatic case, the positive energy part of the Dirac equation for 
the mass eigenstates can be written as 

4f?W* ~ I. (5) 



dt\v 2 J \ E 2 J \v 2 

where = exp(— zcr 2 0) z/ Q (j = 1,2). Integrating Eq.(J5]), we obtain the 
survival probability of the electron neutrino 

P{y e ->■ u e ; R) = cos 2 (o(R) - 9(0)) - sin sin ^ sin 2 ^dr, (6) 

V / 2 2 Jq 2 

where the difference of the eigenvalues is 

AE = E 2 -E X = V(Acos2fl- A) 2 + (Asin2#) 2 . (7) 



In the case of the solar neutrino, since AEdr ^> 1 and |A(0)/A| ^> 1, we 
get 

P(u P -> z/J = sin 2 #. (8) 



2.2 Nonadiabatic transition 

Let us now discuss Eq.(P) in the nonadiabatic case, i.e., in the case where 
we cannot ignore the variation of the effective mixing angle. In this case, we 
have 



dt V u 2 
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The off-diagonal elements in Eq.([n]) stand for an effect in which neutrino 
jumps the energy gap between E\ and E 2 , and such an effect becomes non- 
negligible when the following the adiabatic condition breaks down: 



7 = 



Am 2 



sin 2 29 



2E cos 2$ (ri e /n e ) I 



>1, 



(10) 



which can be derived by comparing the diagonal and off-diagonal elements 
at the resonance point o 

Zener [TT] derived the jumping probability of spins in a linear magnetic 
field. In order for Zener's method to apply for neutrino, the density profile 
of electron has to be approximately linear. In that case, using the adiabatic 
condition Eq. flTU]) . the jumping probability is given by 



Zcncr 



exp 



7T 
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(11a) 



Kuo and Pantaleone [12] showed that the differential equation for the 
exponential profile can be solved exactly. The analytical expression for the 
jumping probability is given by 



exp 


-| 7 (1- tan 2 9) 


— exp 


{-^ 


(\ -tan 2 #y 
V sin 2 # )\ 


1 — exp 


{-b 


fl-tan 2 0V 
V sin 2 9 )_ 





(lib) 



Using the jumping probability P c , we can write down a simple expression for 
the survival probability of the electron neutrino. Again by taking the limits 
f^AEdr > 1 and \A(0)/A\ > 1, we obtain 



P(v P — > v P 



( cos 2 9 sin 2 9 ) 
sin 2 # + P c cos2#. 



Pr 



Pc 
l-Pr 



(12) 



The formulae of Eqs. fllll) for the jumping probability differ when the 
mixing angle is large. In this case, especially for the extreme nonadiabatic 
transition, i.e., 7 — > 0, the jumping probability P c reaches the maximum 
value, i.e., complete conversion, in the Zener's solution, but it does not in 
the exact solution. 



10,. 



1 Here the resonance point differs from the Point of Maximal Violation of Adiabaticity 
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3 Neutrino oscillations with three flavor in 
matter 



The positive energy part of the Dirac equation with three flavor in matter is 
given by 

where a = e, fi, r, the Hamiltonian is 

/ \ / A(t) \ 

H{t) = U A 21 \U 1 + 0, (13b) 
V A 31 / V / 

and the standard parametrization [lj for the 3x3 MNS matrix is 

C12C13 S12C13 Si3e 

-S12C23 - C 12 S 2 3S 13 e l5 C 12 C 2 3 - S 12 S 2 3S 13 e lS S23C13 

S12S23 - c 12 c 23 s 13 e l5 -C12S23 - si 2 c 23 s 13 e lS c 23 c 13 

Here sy = sin^-, = cos 9ij are the mixing angles and 6 is the CP 
phase as in the quark sector [13]. 

3.1 ICA with two nonadiabatic transitions 

For the nonadiabatic transitions with three flavor, the neutrino can have two 
jumping points. When these two points are approximately far apart, from 
the analogy with Eq. (ll2j) . we get the survival probability of the electron 
neutrino HI: 



P{v e -± u e ) = 

= \u el \ 2 p L p H + \u e2 \ 2 p H (i-p L ) \u,,\ 2 a- Pu). (ii) 

where 

1 - Pl Pl 
Pl 1-Pl 
1 

1 
-7/ = 1 1 — Ph Ph 
P H 1-Ph 

Here, Pl and Ph are the jumping probabilities at the lower and higher cross- 
ing point each other, and U e i are the elements of the three flavor MNS matrix. 
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3.2 The setting of our analysis 

What we would like to discuss here is whether the ICA in Eq. (ll4|) holds or 
not. Let us start our study by assuming several situations. 

Throughout this paper, we will take the following reference values for 
the oscillation parameters and the reference function for the electron density 
n e (R) at distance R from the initial point: 



sin 2 26 12 = 0.87 (15a) 

sin 2 2# 23 = 1.0 (15b) 

Am 2 21 = 7.9 x 10~ 5 eV 2 (15c) 

5 = (15d) 

n e (R) ~ 50 x n e (R) cm- 3 (15e) 

n eQ (R) ~245N A exp(-lOR/R & ) cm" 3 , (15f) 



where Na is the Avogadro's number, n e& (R) stands for the electron density 
in the solar standard model [Hj, and R & is the solar radius. 

A few remarks are in order, (i) We refer to the standard value [1] for 
the first three oscillation parameters, (ii) We assume that the CP phase 5 is 
equal to zero. This is because the transition probability of u e — > v e (and its 
antineutrino) did not depend on the CP phase S [T5] . (iii) We assume the 
electron density which is proportional to the one in the Sun but is larger by a 
factor 50, in order to get the two energy level crossing!! For the exponential 
density profile, we can get a simple value n e /n e = const, for Eq. ffTUl) . 

On the other hand, In order to have nonadiabatic transitions at the two 
level crossings, we will have to assume hypothetical values for the following 
parameters: 

sin 2 26 3l 

Am li 
E, 

where we will assume Am 2 x > 0, i.e., we will assume the so-called normal 
hierarchy, throughout this paper for v e (instead of v e ) to have two level 
crossings. 

Furthermore, we adopt the expression Eq.( 111bl) for Pl and Ph for the 
low (high). Here we take each of the adiabatic conditions, 7^ and 7^ as two 

2 If |A/A 3 i| > 1, then it follows that ? e (0) ~ v 3 (0) 



7 



levels at each crossing point. 



From the analogy with Eq. ffTUl) we have 



IL 



A 2 i sin 2 26 



12 



cos2#i 2 {n e /n e ) 1 



(16a) 



A sin 2 26 



13 



cos26>i3 (n e /n e \ 



(16b) 



where we have defined, 



A 2 i = Am 2 1 /2£' 
A31 = Aml 1 /2E 
A = A31 - A 2 i sin 2 6 X2 ■ (17) 

In our study, we use the numerical calculation to check whether ICA holds 
or not. We use the Runge-Kutta method to solve numerically the positive 
energy part of the Dirac equation for three flavor (see Eqs. (lT3]) ). 



3.3 Energy dependence of ICA 

First of all, let us comment on the energy dependence of ICA. Since the 
Hamiltonian in vacuum is inversely proportional to the neutrino energy (see 
Eq. fll3bjl ). the higher the neutrino energy is, the closer the distance of the 
two crossing points becomes. Figure [1] shows the asymptote of the energy di- 
agram, in which the neutrino energy is 50MeV. Although the density profile 
is of the exponential type, the higher point is so near to the lower point that 
the density profile looks like linear. We would like to discuss this problem in 
more detail. 



In FigJ21 we have plotted the probability P{y e —> v e ) by the analytic 
solution assuming ICA and that by the numerical one for each neutrino 
energy. The broken line (the box points) indicates the analytic (numerical) 
solution, respectively. The set of the oscillation parameters assumed here are 

sin 2 26 13 = 0.08 

Am 2 x = 2.4 x 10~ 3 eV 2 . 

We assume the error of the numerical calculation within 0.005, and we find 
no difference between the two solutions in FigfSJ This result indicates that 
ICA holds even if the distance between two crossing points is close, i.e., the 
neutrino energy is high. 
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1 .20585e-1 1 



1.2058e-11 s 



1 .20575e-1 1 



=;1.2057e-11 



1 .20565e-1 1 



1.2056e-11 
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Figure 1: The diagonal elements of the neutrino in matter. Two crossing 
points are close. 
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Figure 2: The difference of P[y e — > u e ) between by the analytic solution 
assuming ICA (the broken line) and by the numerical calculation (the box 
points). Here, sin 2 29 
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From the analogy with Eq.(@J, we get 

(A21 

(A 2 i cos2fl 12 - A) 2 + (A 2 i sin2fl 12 



• 2 n (A 21 sin2fl 12 ) 2 , 1fi v 

sin 20i2 = ti ™ , /a . nn (18a) 



2 % (Asin2fl 13 ) 2 
Sm 2613 ~ (Acos20 13 -A) 2 + (Asin2M 2 ' ( } 

where A has been defined in Eq. (IP7|) . In Fig 131 we have plotted the shapes 
of Eq.f llSaj) — the right line, and of Eq. (I18b|) — the left line, in the case: 

sin 2 26 13 = 1 x 

Am 2 ! = 2.4 x 10" 3 eV 2 
E = 700 MeV, 

as the neutrino propagate in the medium to vacuum. Thus, two resonances 
seem to be far apart. 
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Introducing the notations 



y = A/A 21 
a = A/A 21 . 



we can rewrite Eqs.l TTBj) . 

sin 20i2 = 7 — 7^ , 2 (19a) 



sin 2 2Q X _ 
(cos 2^i2 - y) 2 + sin 2 2#i 2 



9 ~ (a sin20n) 2 . , , 

sin 2 28 13 = - — i — ^1 — -. 19b 

u (acos2# 13 -y) 2 + (asin2# 13 ) 2 v ; 

Thus, each resonance point is 

yn = ot cos 2#i 3 (20a) 
y L = cos2fl 12 , (20b) 

and each half width at half maximums is 

r# = asin2#i3 (21a) 
r L = sin20 12 . (21b) 

These quantities are independent of the neutrino energy. This means that 
two resonances never overlap even though the neutrino energy gets higher, 
i.e., even if the distance between two crossing points gets closer. From this 
one might be tempted to conclude that ICA always holds. The question we 
have to ask here is the validity of the non-overlapping resonance to judge 
whether ICA holds or not. 



3.4 ICA for overlapping resonances 

Let us now imagine the hypothetical situation: 

Am 2 x -»■ 0.04 x Am^ ~ Am 21 , 

that is, 

Am 2 31 /Am 2 21 ~ 0(1). 

This parameter realizes the overlapping resonances (see FigfJ]). Or, let us 
here imagine another hypothetical situation, 

sin 2 2^i3 = 0.75 (0 U = 30°). 
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Figure 5: The difference of P(v e — > v e ) between by the analytic solution 
assuming ICA (the broken line) and by the numerical calculation (the box 
points). Here, sin 2 26 13 = 0.75 and Am^ = 2.4 x l(T 3 eV 2 . 
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This parameter realizes the overlapping resonance too. 

In FigJSl we have plotted the both probabilities P(v e — > v e ) for the second 
situation (sin 2 2^13 = 0.75). Again we find no difference between the two 
solutions. From this we observe that overlapping of the two resonances is 
not a sufficient condition for ICA to break down. 



3.5 Dependence of ICA on the mixing angle and the 
mass squared difference 

To search for the parameter range of ICA, we have examined several cases 
for #13 and Am^, and tested whether ICA held or not. For example, in the 
case: 

Am 2 31 /Am 2 21 ~ 0(1) 
sin 2 2# 13 = 0.08, 

we found that the ICA breaks down (see FigE]). 
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Figure 6: The ICA breaking (Am^/Am^ ~ 0(1), sin 2 2# 13 = 0.08). 
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Table 1: The difference of P[y e — > u e ) between by the analytic solution 
assuming ICA and by the numerical calculation. Nonzero terms mean that 
ICA is broken. 



Am'f.JAmj 2 

sin 2 2<9 13 


1.2 


3.0 


15 


30 


0.75 


0.006 


0.000 


0.000 


0.000 


0.08 


0.048 


0.000 




0.000 


It)" 6 


0.000 






0.000 



In TabJU we have shown the difference between P[y e — > u e ) by the an- 
alytic solution assuming ICA and the one by the numerical one for several 
sets of sin 2 2^13 and Am^/ Am^. The numerical value 0.000 indicates that 
the difference is less than 0.005, the error of the numerical calculations. In 
such cases, we interpreted that ICA holds. The result implies that ICA is 
not applicable when # 13 is large and when Am^/ Am^ is smalljf) 



3.6 The criterion of ICA 



Let us again consider the positive energy part of the Dirac equation for the 
neutrino in medium. From Eq. (113bj) . we obtain 



d _ ~_i ~_ ~ i dll ^ 
dt dt 



(22) 



In the adiabatic case, we can ignore the second term on the right hand side 
of Eq.f l22p . Doing this matrix arithmetic, we have 



i— 

dt 



v 2 {t) 

V h(t) ) 



( 



i0 l2 



\ ic 12 6 



13 



-%9v2 

E 2 



-ic 12 



13 



\ 



^12^13 



E. 



J 



( Mt) \ 

V 2 (t) 



(23) 



3 This is equal to the condition for which the two flavor approximation fails. 
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In order for ICA to apply to Eg. (1231) . the Hamiltonian at the lower 
resonance has to be: 



( E x 

V o 



■i9 l2 \ 
E 2 
£3 / 







(24a) 



Likewise, we have to have the following Hamiltonian Hh at the higher reso- 
nance, for ICA to apply: 



H, 



H 



/ E 1 \ 

E 2 i9 l3 
\ -i0 13 E 3 ) 



(24b) 



In this case, we can get Eq. (1141) . From this, we observe that the ICA breaks 
down, for example, when the extra teams depending on d9\ 3 jdt becomes 
large in Hi at the lower resonance. 

Let us introduce the new parameter for the indicator of the ICA: 



dO 



13 



dt 



AE- 



21 



(25) 



When this parameter is large, it means that the contribution of d9\ 3 /dt is 
large in H L . 

In Figl3 where the x-axis (y-axis) is Am^/ Aro^ (sin 2 2# 13 ), we have 
drawn the contour lines of kl with the neutrino energy satisfying jl = 1 (see 
Eq. liTUl) ). for kl = 1.0, 0.5, 0.1. Thus, the new parameter kl describes the 
results in Tab{T]very well. 



We can rewrite 



k l = A 2 i sin 2^12 / 
A sin 26»i2 



A sin 2 20 13 
A sin 26»i3 



(cos 26-13 - (A21/A) cos26 ) i 2 ) 2 + sin 2 29 13 



cos 29 12 (ri e /n e ) L sin 26*13 

Fixing the neutrino energy by the condition 7^ = 1, we get 



^ (£; = j > = (£) 



[cos 26*13 - (A21/A) cos 26*i 2 ) 2 + sin 2 29 



13 



sin 26*10 sin 26* 



13 



(26a) 



15 




— ' — ' 1 0.001 

A 31 /A 21 



Figure 7: The contour map of K£. 



As in Eq. (|26a|) . we can get k,^ at the higher resonance 

_ Wz? /A 21 \ ((A/A 21 )cos2^ 3 -cos2M 2 + sin 2 2^ 12 
K » {E > lH = 1)= {-^ sin2g 12 sin2fl 13 " (26b) 



Eqs. fl26l) are nothing but the condition of whether ICA holds or not, i.e., the 
condition for ICA as functions of #13 and Am^/Am^. It is remarkable that 
these parameters in Eqs. fl2"6T) are independent of the density profile n e (R). 



4 Disccusion 

In the previous section, we have investigated whether ICA held or not, and 
the neutrino energy which requires ICA to break down ranges in the interval 
5MeV < E < 50TeV. At high energy, the effect of the absorption of the 
neutrino by medium becomes nonnegligible. Let us consider this effect first. 
The reaction number of the neutrino iV(/sec) is 

N = L x <t{uN) 

= N A x p x cr x R 
a(uN) = 0.68 x KT 38 x E (cm 2 ), 

where L is the luminosity of neutrino, a(uN) is the cross section of neutrino- 
nucleus and p is the matter density. Here we estimated the matter density 
p = n e (R) at n e (0) = const. Thus, the flight ranges at each neutrino energy 
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for the present density profile Eq. fll5el) are 

R ~ 2 x KT 1 km 
~ 2 x 10 2 km 
~ 2 x 10 5 km 



(IPeV) 
(ITeV) 
(lGeV). 



Although this is a rough estimate, we find that the neutrino, whose energy is 
not less than ITeV, cannot reach the vacuum area. Therefore, our hypothesis 
that high energy neutrinos are emitted at the center of a supernova and are 
observed outside of the supernova may not make sense. 

Secondly, in dense matter for which there are two crossing points, the neu- 
trino density n u (R) becomes large too. In such a case, the neutrino-neutrino 
interaction is strong HQ H2 US H31 12H1 12H 122] ■ Adding this effect, we need 
to incorporate n v (R), which is time-dependent, into the off-diagonal parts 
for H (see Eqs. (|13p ). No formula of the nonadiabatic transitions attended 
by this effect have been established even the formulae with two flavor. In 
our study, we have assumed that the two jumping points are far apart from 
the range where the neutrino-neutrino interactions become important. 

5 Conclusion 

In this paper, we addressed the question whether the Independent Crossing 
Approximation (ICA) holds or not. Ref. jl] has focused on the picture of the 
non-overlapping resonance by which they judge whether two nonadiabatic 
transitions are independent or not. One of the purposes of this study is to 
check whether this interpretation is correct or not. 

First of all, we have checked the dependence of ICA on the neutrino 
energy. Namely, we checked the dependence of ICA on the distance between 
the higher and lower crossing points. By numerical calculations, we showed 
that ICA is independent of the neutrino energy, that is, the distance between 
two crossing points. 

Secondly, we have checked the validity of the picture of the non-overlapping 
resonance for ICA by varying the parameters #13 and Am^ in the hypothet- 
ical ranges. We found numerically that ICA can hold even if two resonances 
overlap, i.e., ICA does not always break down even if two resonant widths 
overlap. 

Thirdly, we have searched for the case in which ICA breaks down. We 
found that ICA breaks down when # 13 is large and when Am^/Am^ is 
small. 

Finally, we have introduced the new parameters as the criterion of ICA. 
With this criterion of ICA, we interpret the ICA breaking as the sign of a large 
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contribution of the off-diagonal coefficients, for example, the contribution of 
dOis/dt at the lower resonance. We have shown that the new parameters — the 
criterion of 1CA, taking the extra contribution into consideration, describe 
well whether ICA holds or not. 

From the recent experiments, we know that the ratio Am^/Aml is 
large [JJ and the mixing angle #13 is small [23]. Therefore the ICA breaking 
is tiny enough for the two nonadiabatic transitions to be regarded as inde- 
pendent, and all we need is just to multiply Ph and Pl which are obtained 
by the two flavor formalism, as in Ref. [I]. 

In this paper, we have dealt with the only the neutrino-electron interac- 
tion. In other cases, for example, where neutrino transition moments cou- 
ple to the large magnetic fields [251 I2S1 [23 HE], under the influence of the 
neutrino-neutrino interaction [TBI EL3 dSl EES [201 [2D [22], or under the influ- 
ence of the Non-Standard Interaction El [7] , our procedure is necessary. 
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